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Abstract 

We consider a neutral self-interacting massive scalar field defined in a d-dimensional Euclidean 
space. Assuming thermal equilibrium, we discuss the one-loop perturbative renormalization of 
this theory in the presence of rigid boundary surfaces (two parallel hyperplanes), which break 
translational symmetry. In order to identify the singular parts of the one-loop two-point and 
four-point Schwinger functions, we use a combination of dimensional and zeta-function analytic 
regularization procedures. The infinities which occur in both the regularized one-loop two-point 
and four-point Schwinger functions fall into two distinct classes: local divergences that could 
be renormalized with the introduction of the usual bulk counterterms, and surface divergences 
that demand countertems concentrated on the boundaries. We present the detailed form of the 
surface divergences and discuss different strategies that one can assume to solve the problem of 
the surface divergences. We also briefly mention how to overcome the difficulties generated by 



infrared divergences in the case of Neumann- Neumann boundary conditions. 
PACS numbers: 03.70+k,04.62.+v. 
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1 Introduction 



The Casimir effect is the manifestation of the zero-point energy of the quantized electromagnetic 
field, in the presence of metallic plates [1]. A very simple calculation predictes that in a four- 
dimensional spacetime, uncharged perfectly conducting parallel plates should attract each other 
with a force per unit area F(L) oc -p, where L is the distance between the plates. Extensive 
reviews of this subject can be found in Refs. [2] [3] [4] [5] [6]. As stressed by Milloni et al. [7], a 
brief argument showing that the zero point-energy associated with the quantized electromagnetic 
field must have a physical meaning was already given by Einstein and Stern [8]. These authors 
noted that a zero-point energy seems necessary in order to avoid a first-order quantum correction 
to [3~ l in the classical limit (3 » uj in Planck's expression for the average energy of an oscillator 
in equilibrium with radiation at temperature (3~ l . 

Although the vacuum energies of different physical configurations are formally divergent, their 
diference can be finite. In the case of a free scalar field, interacting only with boundary surfaces, 
the Casimir approach can be summarized as follows: first a complete set of modes solutions of the 
Klein-Gordon equation satisfying appropriate boundary conditions, and their respective eigenfre- 
quencies are presented. Next, the divergent zero-point energy is regularized by the introduction of 
an ultraviolet cut-off. Finally, the polar part of the regularized energy is removed using a renor- 
malization procedure. This procedure was first discussed by Fierz [9] a long time ago, followed by 
Boyer [10] and also by Svaiter and Svaiter [11] [12]. In these two last references, an attempt to 
clarify the relation between the cut-off method and analytic regularization procedures in Casimir 
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effect has been developed. In particular, in these papers an analytic regularization procedure was 
interpreted as a cut-off method, and using a mixed cut-off in the regularized zero-point energy, it 
was possible to unify these two methods both in two- and three-dimensional spacetimes. Further, 
a general proof was given that when the introduction of an exponential cut-off yields an analytic 
function with a pole at the origin, then the analytic regularization using the zeta function (or a 
generalization for the zeta function) is equivalent to the aplication of a cut-off with the subtraction 
of the singular part at the origin [13] [14]. More recently, Fulling offered an interesting discus- 
sion with regard to the problems in the renormalization program devised to find the renormalized 
vacuum stress-tensor in different field theories [15]. 

It is important to point out that these results are valid at one-loop level and one is dealing with 
free fields only. It is clear that the formalism must be generalized to take into account the case of 
self-interacting fields. Although higher-loop corrections to the Casimir effect seem beyond experi- 
mental reach today, theoretically such corrections are certainly of interest. Nevertheless, with the 
exception for some few papers, only global issues have been discussed in the study of radiative cor- 
rections to the Casimir effect. One such exception is the discussion presented by Robaschik et al. 
[16]. With this scenario in mind, it is natural to ask the important question: how to implement the 
perturbative renormalization algorithm, assuming the presence of rigid boundaries (hard- walls), 
using the standard weak-coupling perturbative expansion in quantum field theory, that is, how to 
implement the one-loop perturbative renormalization of a self-interacting scalar theory, assuming 
boundary conditions which do break translational symmetry. Our aim when studying these issues 
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is linked to the following question: does the infrared problem have a solution in theories where 
translational invariance is broken? Note that temperature effects can solve the infrared problem 
in some models in quantum field theory [17]; for a recent treatment in non-abelian gauge theo- 
ries at high temperature, and the infrared problem, see for example Ref. [18]. Also, in massless 
scalar \ip 4 theory, if thermal equilibrium with a reservoir is assumed, the infrared problem can be 
solved after a ressumation procedure. The standard is to use the Dyson-Schwinger equation to 
write a non-perturbative version of the self-energy gap equation, or to use the composite operator 
formalism [19] [20] [21]. 

We would like to call the attention of the reader that there are some disagreements in the 
literature as to implementing the one-loop perturbative renormalization in finite size systems 
when translational invariance is broken. In the one-loop approximation, Albuquerque et al. [22] 
found that the mass counterterm depends on the size of the compact dimension in the A</? 4 theory. 
Also, Malbouisson et al. [23] assumed a self-interacting scalar field confined between two infinite 
parallel plates, and using the techniques developed by Ananos et al. [21] these authors didn't 
find any surface countertem in the A</? 4 theory at finite temperature. Furthermore, they were able 
to define temperature and size-dependent mass and coupling constant terms in systems where 
translational invariance is broken. 

The purpose of this paper is to present a detailed calculation of the one-loop renormalization of 
the Ay? 4 theory at finite temperature, assuming that one of the spatial coordinates is confined to a 
finite interval. Since this assumption is not sufficient to explicitly breaking the translational sym- 
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metry, we will further introduce boundary surfaces where the field satisfies appropriate boundary 
conditions. In this situation, the breaking of the translational invariance of the theory is ensured. 
This paper is a natural continuation of the papers of Fosco and Svaiter [24] and also Caicedo 
and Svaiter [25]. Our aim is to further the understanding of the renormalization procedure in 
systems at finite temperature where there is a break of translational symmetry. We will discuss 
the Dirichlet-Dirichlet (DD) and also the Neumann- Neumann (NN) boundary conditions. For the 
Dirichlet-Dirichlet boundary conditions, the model is free of infrared divergences. In the Neumann- 
Neumann boundary conditions case, infrared divergences associated with zero modes will appear 
for bare massless fields. We show that there is no clear meaning for a thermal- or size-dependent 
mass in such situations. Consequently, can not be used to solve the infrared problem in the case 
of Neumann- Neumann boundary conditions a resummation procedure. 

The organization of the paper is the following: in the section II we sketch the general formalism 
of the theory, deriving the one-loop two-point and four-point functions. In section III we use 
two different analytic regularization procedures, i.e, dimensional regularization and zeta-function 
analytic regularization, to identify the polar contributions that appear in the expressions of the 
one-loop two-point and four-point Schwinger functions. In section IV we renormalize the four- 
point Schwinger function and the problem for the infrared divergences is raised. In the conclusions 
we will discuss alternative solutions for the problem of the surfaces counterterms. In this paper 
we use h — c — ks — 1. 
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2 General Formalism and the Finite Temperature Gener- 
ating Functional of Schwinger Functions 

The static properties of finite temperature field theory can be derived from the partition function 
[26]. To obtain the partition function the starting point is the Feynman, Matheus and Salam 
approach [27]. Thus, let us consider the generating functional of (complete) Green's functions for 
a self-interacting scalar field theory defined in a flat d- dimensional Euclidean space Z(h), given by 



where [dip] is a translational invariant measure (formally given by [dip] = l\ xeR ddp(x)) and S[<p] 
is the classical action associated with the scalar field. The quantity Z(h) can be regarded as the 
functional integral representation for the imaginary time evolution operator ((p 2 \U(t2, ti)\ipi), with 
boundary conditions tp(ti,x) = <pi(x) and ip(t 2 ,x) = p 2 (x) which gives the transition amplitude 
from the initial state \<pi > to a final state \ip 2 > in the presence of some scalar source of compact 
support. As usual, the generating functional of the connected correlation functions shall be given 
by W(h) = \nZ(h). In a free scalar theory, Z(h) as well as W(h) can be calculated exactly. 
Regarding the Lagrangian density, we assume that 



where m is the bare mass and Ao is the bare coupling constant of the model. We are also assuming 
rriQ > and also A > 0. The Euclidean n-point correlation functions, i.e., the n-point Schwinger 




(1) 



C(p, dp) = l -{dpf + X -mlp 2 + ^Aoy? 4 , 



(2) 
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functions are given by the expectation value with respect to the weight exp(— S(ip)), defined as 

;,=o- (3) 



G (n) (xiX2 X) = J: ^im. 



The n-point connected correlation functions G^ L \x\,X2- l -,x n ) are given by 

r>M, \ S n W(h) . m 

Finally, the generating functional of connected one-particle irreducible correlation functions (the 
effective action) is introduced by performing a Legendre transformation on W(h), 

r(^ ) = -W(h) + J d d x(p (x)h(x). (5) 

Let us define the proper vertices r^xi, ..,x n ) as: 

r(n) (r r \ _ ^ r (^o) I (a) 

1 [X 1 ,..,X n —- r r w =0, [O) 

where the normalized vacuum expectation value of the field (po(x) is given by 

Mx) = Jh(x)' (7) 

It is clear that in the case of a single scalar field, for a zero normalized vacuum expectation value 
of the field (po(x), the effective action may be represented as a functional power series around the 
value Lpo = 0, with the form 

00 1 r 

r(y?o) = £ -f / A 1 .AjW(x 1) ..,x n )^(xi)...%K). (8) 
„=o n - J 



If the bare coupling constant vanishes, i.e., Ao = 0, the generating functional of all n-point 
Schwinger functions Z{h) can be calculated exactly, since we have to evaluate only Gaussian 
integrals. After some manipulations we obtain that the Gaussian generating functional Z (h) is 
given by 

Z (h) = exp J d d x J d d y h(x) Gg\x - y, m )%)) , (9) 
where the two-point Schwinger function (the inverse kernel) satisfies 

(-A x + m*)G$\x - y, m ) = 5 d (x - y). (10) 

In this situation, the free Euclidean field is a gaussian random variable defined by its two-point 
correlation function 

Gf\x - y, m ) = (x\(-A + mD^y) , (11) 
and the two-point Schwinger function has a well known Fourier representation given by 

0?( -''" J = w/ , W' (12) 

In the next chapter we will show that the two-point function G^\x — y, m ) can be expressed in 
terms of the modified Bessel function of the third kind or Macdonald's function K^(x). At present, 
we are not interested in evaluating the two-point Schwinger function, but only in the analysis of 
the behavior of G ( q\x — y, m ) in a given e -neighborhood. Let us assume that m\x — y\ << 1; in 
this case, for d > 3 we can use that G^\x — y, ml) pz G^\x — y, ml = 0) = \x — y\~( d ~ 2 \ For 
d = 2, the mass parameter can not be eliminated from the denominator and we have the following 
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short distance behavior: G^\x — y,rrify oc ln(m|x — y\). It is well known that a massless two- 
dimensional scalar field theory is not consistent, once the model has severe infrared divergences. 
There are different proposals to circumvent this problem; we only mention some of them. For 
instance, one may violate the positivity of the state vector space; another attempt is to restrict 
the test functions of the theory, and finally one can introduce a cut-off in the definition of the 
positive and negative Wightman functions. It is clear that such cut-off procedure is equivalent to 
introducing a box to regulate the theory in the infrared. Later, we will discusss other strategies 
to solve the problem of the infrared divergences in scalar theories at finite temperature. 

Coming back to the generating functional of all Schwinger functions, for Ao 7^ it is not 
possible to find a closed exact expression for the partition function, and a perturbative expansion 
is mandatory. Let us then assume the weak-coupling perturbative expansion of the theory. It 
is important to point out that the partition function can be defined in arbitrary geometries, and 
classical boundary conditions must be implemented in the two-point Schwinger function, restricting 
the space of functions that appear in the functional integrals. If we want to include thermal effects, 
and assuming thermal equilibrium, from the Feynman, Matheus and Salam formula we have: 

<p b \e- iH ^f-^\<p a ) = r {tf)=iPb exp (1 [ tf dt [ dt^xCfr, dip)) , (13) 



'ip(ti)=(p a \ Jti 

where we have to assume that tf — ij — —if3 and also set (p a = tpb, and the sum over all (p a must 
be performed in order to produce the trace. The partition function Tr e~^ H is given by 

rti-i/3 



Tr 



f [dip]exp (i ( l 1 dt f d d - l xC{^, d(p)) , (14) 

J periodic \ Jti J I 



where the integration over the fields satisfying ip(ti — i(3,x) = <p(t,x). Since the time integration 
must range from some value ti to ti — i/3, let t» = and set the contour along the negative imaginary 
axis from to — if3. Thus, t = —it, where < r < (3, and we have 

Z(h)\ h=0 = [ [d(p]exp( f dr fdt^xCfadip)). (15) 

J periodic WO J J 

To generate the n-point Schwinger functions we need to couples the field with an external source. 
We will assume that the system is confined between two paralel hyperplates, (which we call the 
Casimir configuration), localized at z = and z = L, and we are using cartesian coordinates 
x n _ (r,z), where r is a (d — 1) dimensional vector perpendicular to the z vector. Note that 
since we assume thermal equilibrium with a reservoir, we have periodicity in the first coordinate 
and < r\ < (3. See for example Ref. [28], or for a complete review of quantum field theory 
at thermal equilibrium, see for example Ref. [29]. The choice of Dirichlet-Dirichlet boundary 
conditions means that the scalar field satisfies 

ip(r, z)\ z=0 = (p(r, z)\ z=L , (16) 

and Neumann-Neumann boundary conditions means that 

d d 

-r£P\?, z)\ z=0 = T^Vir, z)\ z=L . (17) 

In the next section we will discuss the perturbative renormalization at the one-loop level of the 
field theory in the presence of rigid boundaries. The great interest of this matter is: when systems 
contain macroscopic structures, how it is possible to implement the renormalization program? 
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We will examinate how does the weak-coupling perturbative expansion and the renormalization 
program can be implemented. In order to identify the singular part of the one-loop two-point 
Schwinger function, we use a combination of dimensional and zeta-function analytic regularization 
procedures. We also present the detailed form of the surface divergences. Note that due to our 
choice (two-parallel hyperplates) , the region outside the boundaries is the union of two-simple 
connected domains. The renormalization of the field theory in such exterior regions must be 
carried out along the same lines as for the interior region. For simplicity we are considering only 
the interior region. 

3 The regularized one-loop two and four-point Schwinger 
functions 

The aim of this section is to reshape a well known result, adding finite temperature effects to 
the problem. In order to implement the renormalization program in a scalar field theory where 
we assume Dirichlet-Dirichlet or Neumann-Neumann boundary conditions on rigid surfaces one 
has to introduce surface counterterms. To write the full renormalized action for the theory with 
rigid boundaries we need two regulators: the first one is the usual e that is introduced in the 
dimensional regularization procedure and the second one which we call r], represents the distance 
to a boundary. Accordingly we will show that the full renormalized action must be given by: 
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+ 



+ 



I d d - x r ( Cl (r?)y2 2 (f, 0) + c 2 ( V )ip 2 (f, L)) 
J d d '\ (c 3 ( V )<p\f, 0) + c A ( V )^(r, L)) 



(18) 



where A(e), B{e) and C(e) are the usual coefficients for the bulk counterterms and the coefficients 
Ci(r)), i = 1, ..4, which depend on the boundary conditions for the field, are the coefficients for 
the surface counterterms. As usual, all of these coefficients must be calculated order by order in 
perturbation theory. Note that we are interested in systems that are invariant under translation 
along directions parallel to the plates, which implies that the full momentum is not conserved. 
For such conditions, a more convenient representation for the n-point Schwinger functions to 
implement the perturbative renormalization is a mixed (p, z) representation. Careless one-loop 
perturbation theory leads to ultraviolet counterterms that depend on the distance between the 
plates or also to the absence of surface counterterms [22] [23]. 

In a straightforward way, in the Matsubara formalism all the Feynman rules are the same as in 
the zero temperature case, except that the momentum-space integrals over the zeroth component 
is replaced by a sum over discrete frequencies. For the case of bosons fields we have to perform 
the replacement 



We begin the study of the interacting theory by building the one-loop correction (G[ ' (A , x, x' 



to the bare two-point Schwinger function G (x,x'), for both the DD and iViV boundary condi- 




(19) 



where we are using the following notation: (/ d d 1 r = J$ dr\ J d d 2 r). 
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tions. Using the Feynman rules we have that G\ (Ao, fi, zi, F 2 , z 2 ) can be written as 

G?\\ , Fi, zi.fb, 22) = y / rf d_1 r dz G 2) (fi - r, *i, z) C? } (0, z) G 2) (f - f 2 , z, z 2 ). (20) 

Even though the functions G^\fi — f 2 , Zi, z 2 ) and G 2 ^(f2 — f 3 , z 2 , z 3 ) are singular at coincident 
points (fi = f 2 , zi = z 2 ) and (f 2 = f 3 , z 2 = z 3 ), the singularities are integrable for points outside 
the plates. It is worth mentioning that the most simple way to take into account the boundary is 
to implement the boundary conditions through the explicit form of the free two-point Schwinger 
function G (x — y, m ). A straightforward substitution yields the order Ao correction to the bare 
two-point Schwinger function in the one-loop approximation for the case of Dirichlet-Dirichlet 
boundary conditions. Using the Feynman rules, G^ 2 \Xq, Xi, x 2 , x 3 , X4) , i.e., the O(Aq) correction 
to the bare one-loop four-point Schwinger function, is given by 

G 2 4 \\o,r 1 ,z 1 ,r 2 ,z 2 ,r 3 ,z 3 ,r 4 ,z 4 ) = y J d d_1 r J d d_ V dz J dz' G 2) (fi - f, z 1: z) 

G ( i\r 2 - f, z 2 , z) (G$\r- P, z, z')f 
G ( *\P - f 3 , z', z 3 )G%\? - f 4 , z', zi). (21) 

Note that we supress the m term in each expression. Again, all Go's are singular at coincident 
points, but the singularities are integrable for points outside the plates, except for G^\f—r', z, z'). 
Having in mind the above discussion, in this section we will study the following expressions: 

^Jd^rJ o L dz (G$\d,z)), (22) 
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and 

^ J d d - l r J d d - l r' £dz J o L dz'(G^\f-r',z,z')) 2 . (23) 

Let us first study \ (0, z) = I(z,m , L, j3,d), and define the following quantities: | = |, 
L = a and finally the dimensionless coupling constant (7 = fi 4 ~ d Xo. Therefore, the argument in 
the integral defined in Eq.(22), I(z, m , a, 6, d) can be written as 

I(z,m ,a,b,d) = — f. 9 V f]sm 2 (—) [d d - 2 p 7 ; 7 . (24) 

There are two points that we would like to stress. First to perform analytic regularizations we have 
to introduce a parameter n with dimension of mass in order to have dimensionless quantities raised 
to a complex power. Second, the generalization for the case of Neumann boundary conditions is 
straightforward, although in this case infrared divergences associated with the n = mode will 
appear in the case of massless scalar field. To circumvent this situation, we must have a finite 
Euclidean volume to regularize the model in the infrared, or trying to implement a resummation 
to generate a thermal mass. We will return to this point latter. 

Using trigonometric identities, it is convenient to write the amputated one-loop two-point 
Schwinger in two parts. The first comprises the contributions that do not depend on the distance 
to the boundary, and the second the contributions that do depend on this distance. Therefore, 
the quantity I(z, m , a, b, d) can be split in two parts Ti(mo, a, b, d) and T 2 (z, mo, a, b, d), i.e.: 

I(z, m , a, b, d) = 21 (m , a, b, d) + T 2 (z, m , a, b, d). (25) 

The first quantity Ti(m , a, b, d), independent on the distance to the boundaries can be expressed 
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in the following way: 

Ti(m , a, 6, d) = Io(m , a, b, d) + h(m , a, b, d) + 7 2 (^0j a , d), 
where each term is given respectively by: 

Jo(m 0) a, 6, d) = - 16(2 ^_ 2q6 / ^(^g) » 

g 00 f 1 

A (m , a, 6, d) = d _ 2 E / ^Pj^Z 2 , , BM2 v 

8(2tt) a6 n=1 ^ (p 2 + m 2 , + (^) 2 J 

and finally 

The contribution that depends on the distance to the boundaries given by T 2 (z,m , a,b, 
be split in the following way: 

T 2 (z, m Q , a, b, d) = I 3 (z, m , b, d) + I A (z, m , a, b, d) + I 5 (z, m , b, d) + I 6 (z, m , a, b, d). 

Each term contributing to T 2 (z, m , a, b, d) is given, respectively by: 

I 3 (z, m , b, d) = -zrh(d) / dv(v 2 — m )~ exp(— 2v z), 

2b J mo 

l4(z,m , a,b, d) = —h(d) / dv (v 2 — ml)~ (cothav — 1) cosh2vz, 

2b J mo 

I 5 (z,m ,b,d) = -h(d)^2 dv W - m l - ("r) 2 exp(-2vz), 

n=1 Jmo V / 

and finally 

I 6 (z,m ,a,b,d) — -h(d)^2 / dv I v 2 — ml — (— ) 2 ) (cothaw — 1) cosh2vz. 
b n=1 J a \ b J 
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In the above expression the quantity a is given by 



a=(ml + (^ry\ (35) 

and h(d), that appears in Eqs.(31), (32), (33) and (34) is an entire function given by 

Kd) = ^ ^ (36) 
4(4tt) 2 T(— ) 

Let us investigate each contribution in detail. Using dimensional regularization we obtain for 
Io(mo,d) the following expression: 



An analytic expression for the Gamma function T(z), defined in the whole complex plane, can 
be found and in the neighborhood of a pole z = —n, (n — 0, 1, 2..) the Gamma function has the 
representation 

(-l) n 1 

T(z) = ^-^- + n(z + n), (38) 

n\ (z + n) 

with regular part Q(z + n). Using that 4 — d — e an the duplication formula for the Gamma 
function F(z) we have 

I (m ,a,b,d)\ d=4 = (- + n(e)) . (39) 

Here one may adopt different renormalization schemes. We can choose the minimal subtraction 
(MS) scheme, in which we eliminate only the pole term ^ in the dimensionaly regularized expression 
for the Schwinger functions. Another choice is the modified minimal subtraction (MS) scheme, 
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where we eliminate not only the pole term i but also the regular part around the pole. Note that 
in the minimal subtraction scheme the counterterms acquire the simplest expression, while the 
renormalized Schwinger functions have more complicated expressions. Let us analyse the second 
expression, given by h(mo, a, b, d). Using dimensional regularization it is possible to show that 

a d 00 1 

h ^ °' 6 ' $ = R(9 A*-2 nh T ( 2 " ?) £ ~ ( 4 °) 

8{2y/n) d 2 ab 2 „=i L2 + (»)2 ^ 

We note that to extract a finite result from h(m , a, b, d) we still have to use the analytic extension 
of the Epstein-Hurwitz zeta function. A direct calculation gives 

^ d s 



Mmo ,a, M ) = -_^„ i « I _^ rT r ( 2--) + 

+ 4 (am n n) ^ K 3 -d (2m n na) ) . (41) 

n=l 



gm$- 3 1 / /3-d N 



86 {2n) d - 1 \ \ 2 1 — 



The first term in the above equation is a polar part and the second one is finite. Assuming the 
minimal subtraction scheme, h(mo, a, b, d) becomes finite. The next term that we have to analyse 
is hiirio, a, b, d) defined by: 

a 1 00 r 1 



I 2 (m ,a,b,d) = -f- . . . - Y d d - 2 Pl r-. (42) 

The contribution given by the above equation is a part of the amputated one-loop two-point 
Schwinger function that does not depend on the distance to the boundaries, and in the renormal- 
ization procedure it will require only a usual bulk counterterm. The form of the counterterm is 
given by the principal part of the Laurent expansion of Eq.(42) around some d, which must be 
given by the analytic extension of the Epstein zeta function in the complex d plane. The structure 
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of the divergences of the Epstein zeta function is well known in the literature [30] [31] [32] [33]. 
Since the polar structure of the above equation can be found in the literature, we will focus only 
on the position-dependent divergent part given by T 2 (z,mo,a,b,d). We are now in position to 
discuss the behavior of I${z, m , b, d), h(z, m , a, b, d), h(z, m , b, d) and finally Iq(z, mo, a, b, d). 

Let us first analyse I 3 (z,m , b, d). Using the following integral representation of the modified 
Bessel functions of third kind, or Macdonald's functions K v {x) [34], 

r°° 1 9ti 1 

/ (x 2 -u 2 y- 1 e - fix dx = -==(— ) v -^r(u)K v _i( Uf i), (43) 

Ju y 7T /i 2 

which is valid for u > 0, Re (fx) > and Re (u) > 0, we see that Is(z, m , a, b, d) can be written in 
terms of these functions. A simple substitution gives 

r t , ,n 2 hid) .mn.d-3 T ^ . . . 1 1 , 

I z (z, mo, a, b,d) = - ^/=p=i (— ) 2 K^(2m z). (44) 

Using a asymptoptic formula for the Bessel function, I${z, mo, a, b, d) is given by 

/ 3 ( z ,m ,a, M ) = g p^Effi (45) 

We can see that we have a divergent behavior as z — > 0, which demands a surface counterterm. 
Let us show that the other terms also contain surface divergences, and study h(z, mo, a, b, d). To 
advance in the calculations, we have to extend the binomial series for both positive or negative 
integral exponents, written in the form 

oo 

(l + x) k = J2C k n x n . (46) 

n=0 
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First, it is possible to show that the binomial expansion holds for any real exponent a, \x\ < 1 
and ae R, i.e., 

oo 

(l + x) a = '£CZx», (47) 

n=0 

where C™ are the generalization of the binomial coefficients. Since we are using dimensional 
regularization, it is possible to extend the binomial expansion when both the exponent a as well 
the variable x assume complex values. For this purpose we use the following theorem: 
For any complex exponent a and any complex z in \z\ < 1, the binomial series 

oo 

J2C:z n = l + C l a z + .. + C:z n + .. (48) 

n=0 

converges and has for sum the principal value of the power (1 + z) a , where the principal value of 
the power b a is given by the number uniquely defined by the formula b a = exp(alnfr), where a and 
b denotes any complex numbers, with b ^ as the only condition, and In b is given its principal 
value. Going back to J 4 (,z,m , a,b, d), using the generalization of the binomial theorem, let us 
define C^(d,k) = Ih^-lfC^ to obtain 



h(z,m ,a,b,d) = -^7j2 cW ( d ' k )( m o°>) 2k / « d " 4 " 2fc (coth W -l)cosh(— ). (49) 

Let us use the following integral representation of the Gamma function, 

f°° I 

/ dtt^e-^ = —T(fi), Re{v)>0, Re(u) > 0, (50) 
jo W 

and also the following integral representation of the product of the Gamma function times the 
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Hurwitz zeta function 

/ dt^- 1 e _at (cotht-l) = 2 1 -' 4 r(/i)C(/i,7r + l) i2e(a) > 0, ite(/x) > 1, (51) 
Jo 2 

where ((s, u) is the Hurwitz zeta function defined by [34] 

oo -y 

C(g,if) = E ? j_ w #e(s)>l 0,-1, -2... (52) 

n=0 V 71 + W / 

It is not difficult to show that I^z, m , a, b, d) contains surface divergences at z = and also 
z = a. For more details, see for example Ref . [35] . The other expression that we have to study is 
Is(z, m , a, b, d). Using an integral representation of the Bessel function of third kind we have: 

I 5 (z,m ,a,b,d) = \ — j_ ^(-^Kd-s (2az). (53) 
b (200* n=1 z 

Using an asymptotic representation of the Bessel function it is posssible to present also the singular 
behavior near z — 0. Let us finally investigate I 6 (z,m ,a,b,d). A simple calculation for the 
massless case gives 

I 6 (z, m , a, b, d)L=o = -5Z37 £ C (2) (4 k)(-) 2k £ ^ / d« u d ~ 4 - 2k (cothu - 1) cosh( ), 

CL j r> „ 1 «/ — ; — 

/c=0 71=1 b 

(54) 

where C^(d,k) = h(d)(—l) k C k - 4 n 2k is an entire function in the complex d plane. The integral 

2 

that appear in Eq.(54) cannot be evaluated explicity in terms of well known functions. Nevertheless 
it is possible to write Eq. (54) in a convenient way where the structure of the divergences near the 
plate when y — > b appear. Clearly for details see Ref. [35]. In the next section we will investigate 
the singularities of the four-point Schwinger function. 
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4 The four-point Schwinger function in the one-loop ap- 
proximation 

We now turn our attention to the four-point Schwinger function in the one-loop approximation. 
For simplicity we shall study only the zero temperature case. In this section we are following the 
discussion developed in Ref. [25]. Introducing new variables as u± = z ± z', and also (p — f — r'), 
the zero-temperature two-point Schwinger function in the tree-level G (p, z, z') can be split into 

Gf (p, z, z') = G?(p, u + ) + G m (p, u_), (55) 

where we are defining A n (a, mo, d, p) by 

1 r 

A n (a,m ,d^ = j-^J d^p {p + {rf)2 + miy (56) 

and so G±\p, u±) can be expressed as 

G ( ±{p,u±) = =p- cos(^^) A n (a,m ,d,p). (57) 

Before proceeding, let us present a explicit formula for the free two-point Schwinger function 
G±\p,u±) in terms of Bessel functions. Let us define an analytic function f(d) by 

1 r(4=i) 

M = -=71 ^ r /a • (58) 
v/7r(27r) 2 1 {—) 

(2) 

Strictly speaking, it is possible to show that we can write G± (p, u±) in terms of the Bessel function 

of third kind. To this end, we use the standard formula 

1 r - 1 rfstl^ roc j 

— — / d d rF(r)e ik - r = / F(r)r2 J d - S (kr)dr, (59) 

(2n) d J v^(27r)f r(^)7o — V ; K J 
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which leads us to: 

G±'{P,u±) = T^r- ^ cos(— — ) / dp —————— J^(pp), (60) 

p 2 a n= i a Jo [p + {—) + m ) 2 

where J^x) is the Bessel function of the first kind of order v. The integral in Eq.(60) can be 
calculated by using the result [34] 

f t 9w^ lf ' w ' (61) 

implying that it is possible to write G±\p,u±) as 

CL 2 '(P, u ± ) = £ cos(=±) ((™)» + mi)*K v ( P f^?) ■ (62) 

p 2 a n =i a \ a / 2 \ v a/ 

Using Eq.(55) and the above formula, the explicit expression for the two-point Schwinger function 
in a generic (/-dimensional Euclidean space confined between two flat paralel hyperplanes, where we 

assume Dirichlet-Dirichlet boundary conditions is given. It is difficult to use the above expressions 

(2) 

for G± (p,u±) to investigate the analytic structure of the four-point function for both the bulk 
and near the boundaries. Nevertheless, it is clear that the divergences of the four-point function 
in the one-loop approximation appear at coincident points and therefore the singular behavior is 
encoded in the polar part of M(X , a, m, d) given by 

M{X ,a,m ,d) = g 2 J d^r J d d ~ x r' £ dz £ dz' F{r,r< , z, z') (G { 2) (f - r> , z, z')f . (63) 

It is easy to show that G^\\o, a, m , d) amp is given by 

G 2 (Xo,a,m ,d) 

amp 

2 r r r r 00 pip-(q-k) 

, o 9 / d d -\ \ d d -\' / d d ~ x k / d d -\y^ — =, , (64) 

2(27r) M - 2 i J J J H ti (^ + (f) 2 + mg)(F + (f) 2 + mg)' V 
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where F(f, r r , z, z') is a regular function. As with the one-loop two-point function, it is not difficult 
to realize that the above equation has two kinds of singularities, those coming from the bulk and 
those arising from the behavior near the surface. As before, the behavior in the bulk is similar to 
the thermal field theory case and consequently we will discuss only the singularities arising from 
the boundaries. This can be done studying the polar part of M(X Q , a, mo, d) given by 

M(\,a,m ,d) = ?- dz dz'T{z, z') (g^ (0, z, z')) , (65) 

Zi J J 

where J-(z, z') is a regular function. Now, we recall that the form of G±\p, u±)\ p=0 is given by, 

Ofi^^^t^j^ ^J^^ (66) 
where it is not difficult to show that 

G±\p,u±)\ p=0 = =F ^-^-A ( y p,L,m )\ p= o + f2(a,m ,d, . (67) 
In the above definition we are making use of the auxiliary function f2(a, d, m , z) given by 

f 2 (a,m ,d,z)= d V-f—= 2 . ~ r~2~ji 2\ • (68) 

J Jp 2 + mo smh(aJp 2 + mg) 

Note that the amputated one-loop two-point Schwinger function can be decomposed in a transla- 

tional invariant part and a translational invariance breaking part, given exactly by /^(a, mo, d, z). 

When we sum to find the free propagator, we end up with the following expression 

Go\p,z,z')\ p=0 = f 2 (a,m ,d,^=-) - f 2 (a, m , d, (69) 
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For the sake of simplicity, we will discuss only the massless case once the singularities of the massive 
case have the same structure as in the massless one. The function /2(a, m , d, is non-singular 
in the bulk, i.e., in the interior of the interval [0,a], while / 2 (a, m ,c?, ^) has a singularity along 
the line z = z' . Indeed, closer inspection shows that for < z, z' < a the only singularities are 
those at u + = 0, u+ = 2a and also U- = 0. The former two are genuinely boundary singularities 
(the two conditions imply z, z' — > or z, z' — > a), while the last comes from z = z' in the whole 
domain and is just the standard bulk singularity. In fact, using the structure of the two-point 
function and showing just those terms from which singularities might arise, one finds that the 
counterterms for M are given by 

-pole j Q dz I dz'[ I - T ± ¥ - 2 + - -1 + ^— JL— + (70) 

where Q,i = 1,..3 are regular functions that do not depend on z or z' . From this discussion 
it is clear that in order to render the field theory finite, we must introduce surface terms in the 
action. This is a general statement. For any fields that satisfy boundary conditions that break the 
translational invariance it is suffices to introduce surface counterterms in the action, in addition 
to the usual bulk counterterms, to render the theory finite in the ultraviolet [36] [37] [38]. Now 
we are able to discuss whether in the Casimir configuration the infrared problems can be solved 
for the case of Neumann boundary conditions. For the case of massless (X^> 4 )d theory at finite 
temperature, the infrared problem can be solved after a resummation procedure [17] [18] [19] [20] 
[39]. The key point for the solution of the infrared problem is to use the Dyson- Schwinger equation 
to rewrite the self-energy gap equation. Simple inspection of Eq.(24) show us that it is not possible 
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to implement such scheme in a situation where there is a break of translational invariance. 

A different possibility to approach the infrared problem is to single out the zero mode compo- 
nent of the field, treating the non-zero modes perturbativelly and treating the zero mode exactly. 
This is a standard procedure in high-temperature field theory, where by means of the dimensional 
reduction idea, we relate the thermal Schwinger functions in a d-dimensional Euclidean space to 
zero temperature Schwinger functions in a (d — 1) dimensional Euclidean space [40] [41] [42]. In 
this situation we have a dimensionally reduced effective theory. The key point in this construc- 
tion is the fact that the leading infrared behavior of any field theory at high temperature in a 
d-dimensional Euclidean space is governed by the zero frequency Matsubara mode. 

5 Discussions and conclusions 

In this paper we were interested in the analysis of the important questions of perturbative expan- 
sion and renormalization program in quantum field theory with boundary conditions that break 
translation symmetry, assuming that the system is in equilibrium with a reservoir at tempera- 
ture Specifically, the purpose of this paper is to study the renormalization procedure up to 
one-loop level in the (A^ 4 )^ theory at finite temperature assuming that the scalar field satisfies 
Dirichlet-Dirichlet or Neumann-Neumann boundary conditions on two parallel hyperplates. 

We first obtained the regularized one-loop diagrams associated with scalar field defined in 
the Casimir configuration in a d-dimensional Euclidean space. We obtained a well-know result 
concerning surface divergences that appear in the one-loop two-point and four-point Schwinger 
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functions as a consequence of the uncertaintly principle. There are at least three different possible 
solutions to eliminate these divergences. The first one is to take into account that real materials 
have imperfect conductivity at high frequencies. As was stressed by many authors, the infinities 
that appear in renormalized values of local observables for the ideal conductor (or perfect mirror) 
represent a breakdown of the perfect-conductor approximation. A wavelength cutoff correspond- 
ing to the finite plasma frequency must be included. The second one is is to substitute classical 
boundary conditions by classical potentials; for previous papers using this idea see for example 
[43] [44] [45]. A localized boundary with some cut-off can also be used to replace the potential. 
Nevertheless, it is necessary to renormalize the potential [25]. The third one regards a quantum 
mechanical treatment of the boundary conditions . A fruitful approach to avoid surface diver- 
gences, discussed by Kennedy et al. [46] is to treat the boundary as a quantum mechanical object. 
This approach was developed by Ford and Svaiter [47] to produce finite values for the renormalized 
< ip 2 > and other quantities that diverge as one approach the classical boundary. 

Consequently, we have two main distinct directions for future investigations. The first is related 
to the infrared divergences of our model. Infrared divergences of massless thermal field theory arise 
from the zero frequency Matsubara modes, so we construct an effective (d— 1) dimensional theory 
by integrating out the nonstatic modes and therefore the zero frequency Matsubara modes which 
are responsible for infrared divergences can be treated separately. The second direction is related 
to the surface divergences. In the Euclidean formalism for field theory, one may imagine that 
our simplified model of rigid boundaries is a good approximation only for points in the bulk; for 
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points close to the surfaces however, our approximation is no longer acurate and a model taking 
into account at least thermal fluctuations of the boundaries must be developed [48]. In other words, 
a fundamental understanding of the perturbative renormalization algorith in the standard weak- 
coupling perturbative expansion of an Euclidean field in the presence of fluctuating boundaries is 
desired. This interesting situation of thermal fluctuating boundaries is under the investigation by 
the authors. 
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